1 



o 
o 

c 

in 

m 

< 

u 
S 



in 

> 

On 
O 
(N 

O 

-i— > 

£ 

> 

X 



g-Classical orthogonal polynomials: A general difference calculus approach 



R. S. Costas-Santos (rscosa@gmail.com) 

Department of Mathematics. University of California, Santa Barbara, CA 93106, USA. 
F. Marcellan (pacomarc@ing.uc3m.es) 

Departamento de Matemdticas. Universidad Carlos III de Madrid. Avenida Universidad 30, 28911-Leganes, 
Spain 

June 23, 2009 

Abstract. It is well known that the classical families of orthogonal polynomials are characterized as the 
polynomial eigenfunctions of a second order homogenous linear differential/difference hypergeometric operator 
with polynomial coefficients. 

In this paper we present a study of the classical orthogonal polynomials sequences, in short classical OPS, in a 

more general framework by using the differential (or difference) calculus and Operator Theory. 

The Hahn's Theorem and a characterization Theorem for the g-polynomials which belongs to the g-Askey and 

Hahn tableaux are proved. Finally, we illustrate our results applying them to some known families of orthogonal 

^-polynomials. 
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1. Introduction 

Classical orthogonal polynomials constitute a very important and interesting family of special 
functions. They are mathematical objects which have attracted attention not only because of 
their mathematical value but also because of their connections with physical problems. In fact, 
they are also related, among others, to continued fractions, Eulerian series, elliptic functions 
[7, 13], and quantum algebras [17, 18, 28]. 
They also satisfy a three-term recurrence relation (TTRR) [24] 

x{s)p n (x(s)) = a n p n+1 (x(s)) + f3 n Pn{x{s)) + 7„p n _i(x(s)), n > 0, 

where 7„ / 0, n > 1, as well as their derivatives (or differences or (/-differences) also constitute 
a a sequence of orthogonal family (see e.g. [3, 4, 11, 24] for a more recent review). 

Indeed, a fundamental role is played by the so-called characterization Theorems, i.e. the 
Theorems which collect those properties that completely define and characterize the classical 
orthogonal polynomials. 

One of the many ways to characterize a family (p n ) of classical polynomials (Hermite, Laguerre, 
Jacobi, and Bessel), which was first posed by R. Askey and proved by W. A. Al-Salam and T. 
S. Chihara [1] (see also [21]), is the structure relation 

(p(x)p' n (x) = a n p n+ i{x) + l n p n {x) + C n p n -i(x), 71 > 0, (1) 

where <p is a fixed polynomial of degree at most 2 and c n ^ 0, n > 1. 

A. G. Garcia, F. Marcellan, and L. Salto [14] proved that the relation (1) also character- 
izes the discrete classical orthogonal polynomials (Hahn, Krawtchouk, Meixner, and Charlier 
polynomials) when the derivative is replaced by the forward difference operator A defined as 

A/(*) := {£+ - S)f(x) = f(x + 1) - f(x), 



where S + is the forward shift operator and J? the identity operator. 
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Later on, J. C. Medem, R. Alvarez-Nodarse, and F. Marcellan [22] (see also [2, 4, 5]) char- 
acterized the orthogonal polynomials which belong to the g-Hahn class by a structure relation 
obtained from (1) replacing the derivative by the g-difference operator <2> q defined as 

The aim of the present paper is to go further in the study of classical polynomials giving 
two new structure relations for the g-polynomials which belong to the g-Askey tableau and 
the lattice is g-quadratic, i.e. x(s) = c\q s + C2q~ s + C3, with c\C2 7^ 0, being q £ C := C \ 

({0}u(U n >i{«€C:^ = l})). 

In fact, we prove that the following relation characterizes the g-polynomials. 

yy ( , « A Pn+l(g) . , Ap TO (s) Ap ra _l(fi) 

^(x(s)) = e n ^^ + f n ^- +9n ^^-, n>0, 
where ^# is the forward arithmetic mean operator: 

and, (e n ), (f n ), and (g n ) are sequences of complex numbers such that e n 7^ and g n 7^ 7„. 
Furthermore, we give a characterization Theorem for classical orthogonal polynomials in a more 
general framework using Operator Theory, as well as we deduce some unified expressions for the 
second order homogeneous linear differential (or difference) hypergeometric operator and for its 
polynomial eigenfunctions. 

The structure of this paper is as follows: In Section 2 we introduce some notations and definitions 
used throughout the paper. In Section 3 we define the Rodrigues Operator and using Operator 
Theory we deduce a unified expression for the linear differential (difference resp.) hypergeometric 
operators related to the classical families, and for their polynomial eigenfunctions, as well as 
other algebraic properties. In Section 4 we present a new characterization Theorem for g-classical 
orthogonal polynomials and we prove the Hahn's Theorem for the g-polynomials of the g-Askey 
tableau. In Section 5 we illustrate our results applying them to some known (/-polynomials. 



2. Preliminaries 

The standard classical orthogonal polynomials (Hermite, Laguerre, Jacobi, and Bessel) are 
eigenfunctions of a second order linear differential operator [26] 

^:=^)^+f(x)^, (2) 

where a and f are polynomials of degree at most 2 and 1, respectively. This operator is also said 

to be a standard Hamiltonian operator (see [24]). 

In fact, a sequence of monic polynomials, (p n ), such that 

Jfp n {x) = (^nf' + n(n- l)^jp n (x), n > 0, 

is orthogonal with respect to a certain weight function p{x) supported on Q C R, i.e. 

Pn(x)p m (x)p(x)dx = dl5 ntm , n, m > 0. (3) 

Jn 

Notice that £1 is an interval associated with the regular singularities of Jtif. 
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Moreover, this weight function p(x) fulfills the following Pearson equation 

^-(a(x)p(x)) = f(x)p(x). 

On the other hand, we can consider two different discretizations of the operator (2) replacing the 
derivative by certain approximations on a lattice. Thus we need first to introduce the concept 
of lattice. 

DEFINITION 2.1. [11, 24] A lattice is a complex function x G C 2 (A) where A is a complex 
domain, No C A, and x(s), s = 0, 1, . . . are the points where we will discretize (2). 

For a first approximation of the operator , we consider the uniform lattice x(s) = s. Thus the 
following second order homogeneous linear difference hypergeometric operator 

:= cr(s)AV + t(s)A, (4) 

can be introduced where V and A are the backward and forward difference operators, respec- 
tively, being 

V/(x) = (g- - jr)f[x) = f(x) - f{x - 1), 

a(x) := a(x) — \f{x) is a polynomial of degree at most 2, and r(x) = f(x). 
Indeed, a sequence of polynomials (p n ) satisfying 

J^APn(x) = (nf' + n(n - l)-y) Pn(x), n > 0, 
is orthogonal with respect to a certain weight function p{x) supported on SI C I, i.e. 

^2 Pn{x)p m {x)p{x) = d 2 n 5 n , m , n, m > 0, (5) 

with some extra boundary conditions (see §2.3 in page 26 of [24]). 

As in the continuous case, the weight function p fulfills the following A-Pearson equation 

A(a(x)p(x)) = r{x)p(x). 

The second way to discretize the operator is to consider a nonuniform lattice x(s), (see [24]) 
which leads to the following linear difference hypergeometric operator: 

A V A 

Jt? g := <j(s)— — r^ =—-r + r ( s H — (6) 
q v 'Vxi{s) Vx(s) v ' Ax(s) v ; 

where x(s) = c\q s + C2q~ s + C3, q G C, a(s) := a(x(s)) — |r(x(s))Vxi(s), x±(s) := x(s + \), and 
t(s) = f(x(s)). 

Notice that, in general, o~(s) is not a polynomial on x(s) but if either c\ = or C2 = then 
a(s) is a polynomial on x(s) of degree at most 2, i.e. a(s) = a(x(s)). 
Indeed, a sequence of monic polynomials (p n ) satisfying 

J^ q p n (x(s)) = [n] q (a q {n - 1) f' + [n - l] 9 yj p n (x(s)), 11 > 0, 

where [s] q denotes the q- number (in its symmetric form) 

qi — q~z 



s e C, 

q2 — q 2 
and 



<M S ) := 2 ' S ' 



is orthogonal with respect to a certain linear functional u G P', i.e. 
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(u,p n p m ) = d 2 n 8 m:n , m,n>0. 



REMARK 2.1. Notice that u admits a Lebesgue-Stieltjes integral representation, not necessarily 
unique, 

<u,p) = J^p(z)dp(z), 

where T is a contour in the complex plane and p, is a measure, with supp(p) = T. 

In fact, in most of the cases, such an integral can be transformed into a sum in such a way 
the linear functional can be written as [16] 

(u,p) = P( x ( s i))P( s iW X l( s i)i S i+l = S i + !) ( 7 ) 

where the function p is supported on f2 C BL For this reason we will assume the above represen- 
tation for the functional although the proofs when u can not be represented in such a way are 
analogous and we will omit them. 

Moreover, as in the discrete case, the weight function p(s) satisfies the following (/-Pearson 
equation 

A(a(s)p(s)) 



Vxi(s) 



r(s)p(s). (8) 



REMARK 2.2. It is important to point out that in the orthogonality conditions (3), (5), and 
(7) one needs to impose some boundary conditions, which follow from the fact that the moments 
of the corresponding measure must be finite (see, for instance, [12], [24, p. 27, p. 70], [25]). 

An important fact related to these operators is its hypergeometric character, i.e. if y is an 
eigenfunction of the linear operator J>$? , J%a, or Jt? q respectively, its derivative y = y', its 
difference y = Ay, or its (/-difference y = A^y, respectively, is an eigenfunction of the second 
order linear operator of the same kind. For instance, y = A^y is an eigenfunction of the 
following second order linear difference operator: 

■= -(^)^TT^T + U + 1)^# + ^ + 1} " a{s) 



Vxi(s)Vx(s) V Ax(s) Ax(s) / Ai(s)' 

where 

A (m) := . A . r . A -. , ...-r^r, m>l, x m {s):=x(s + f), m G Z. 
Ax m _i(s) Ax m - 2 (s) Ax(s) 1 

Throughout the paper we will study some properties related to the difference operator and 
its eigenfunctions. 

We can rewrite ^ as follows 

Moreover, with this definition we can also rewrite the (/-Pearson equation (8) as 

a(s + l)p(s + 1) = (a(s) + r(a)Vxi(a))p(a). (10) 
Combining (9) and (10), we get the symmetric or self-adjoint form of (6) 



1 V <\ 



p(s) Vxi(s) 



A 

o 



A X ( s y 



(ii) 



where po(^) := p(s) and pi(s) := p(s + l)a(s + 1). 

This representation for the difference operator (6) is the key to define the Rodrigues Operator 
which we will use in order to unify all the representations in the context of the classical orthogonal 
polynomials. 
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REMARK 2.3. A direct calculation yields the analogous formulas for the operators J4? andJt?&. 



1 d 



_p(x) dx 

where po(x) := p{x) and p±(x) := p{x)a{x); and 



Pi 0*0 



P(s) 



Vpi(s) 



d_ 

dx ' 



oA, 



where po(s) := p(s) and pi(s) := p(s + l)er(s + 1). 

Using this information we can write the polynomial eigenfunctions of the operator Jif q as (see 
[24, (3.2.10), p. 66]) 



p n (x(s)) :-- 



Br, 



Pi( s ) 



_po(s) Vxi(s) 
or, after some simplifications, as 



V 



Pi(s) Vx 2 (s) 



P2(s) 



V 



Pn-l(s) Vl n (s 



Pn{s) 



, (12) 



Pn(x(s)) 



B 



n ^(n) 



Pols J 



where po(s) = p(s), Pk(s) = pk-i(s + l)cr(s + 1), for every k £ Z, and 



v w = v w v 
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V 



V 



V 



Va; fc+ i (s) Vx k+ i (s) Vx n (s) ' 



fe = 0, 1, . . . , n — 1. 



2.1. The (forward) mean arithmetic process 

It is well known that for any pair of polynomials, ir and £, the following relations hold: 

A^faO = (A^tt^ + 

A (1) «) = (A«£H + (^ + 0A (1 V. 

We will say that we apply the (forward) mean arithmetic process to AW(7r^) if we get the 
arithmetic mean of the above two expressions obtaining in such 

A^K) = (A«tt)^ + (A«£)^#tt. 

For example, 

A^(x(s)p n (x(s))) = J? Pn (x{s)) + aq {l) Xl {s)^p n {x{s)) + c 3 (q)(l - a q (l)). 

In an analogous way we define the (backward) mean arithmetic mean process although we will 
not use it. 



3. The Rodrigues Operator 

The expressions given in the previous section for the difference operator Jf? q and for its poly- 
nomial eigenfunctions (12) suggest us to consider a new operator which we call the Rodrigues 
Operator. 

DEFINITION 3.1. Given functions a and p, where p is supported on Vt, and a lattice x(s), we 
define the k-th Rodrigues Operator associated with (a(s), p m (s), x m (s)) as 

^o(o-(s),pm(s),x m (s)) := y, Mi(o(s),p m (s),x rn (s)) := — V p m+ \{s) o/, 

Pm( s ) v X m Jf-l{S ) 

M k (a(s),p m (s),x m (s)) := &i(o-(s),p m (s),x m (s)) o^fc_i(<T(s),p m+ i(s),x m+ i(s)), k > 2. 



For the sake of simplicity, we denote 

@k{Pm(s),x m (s)) :=@ k (a(s),p m (s),x m (s)), m G Z, fe > 0. 
Notice that in the continuous case 

@l(pm(x),x) = l —-^-p m+1 (x) o J, 

p m (x) dX 

and in the discrete case 

&i(Pm(s), x(s)) = — l —-Vp m+ i{s) o J. 

Pm\ s ) 

From this definition the difference operator (11) can be written as 

^ = ^i(p(s),s(s))oAW (13) 

and (12) reads 

Pn {x(s)) = B n M n {p{s),x(s)){l), B n ^0, n>0. 

Therefore, this is the way to write the operators Jt? (J^a or J4? q , resp.) and their corresponding 
polynomial eigenfunctions in a unified way. 

Now, let us consider other properties related to the Rodrigues Operator. 
LEMMA 3.1. The Pearson-type equation (8) can be written as follows 

<% 1 (o-(s),p(s),x(s))(l)=T(s). 

As a consequence we get result: 

LEMMA 3.2. [24, (3.2.18), p. 66] For every integers n, k, < k < n, there exists a constant, 
C n! k, such that 

A^p n (x(s)) = C nik <% n _ k (p k (s),x k (s)){l). 



4. The characterization Theorem for the ^-classical polynomials 

The characterization Theorems constitute a very useful tool to analyze classical orthogonal 
polynomials as we stated in the introduction. Recently, in [4] and [2] a characterization of q- 
classical polynomials in the ^-linear lattice x(s) = cq ±s + d is presented. We will extend these 
results for the (/-quadratic lattice 

x(s) = ciq s + c 2 q~ s + c 3 , q G C. (14) 

For our purpose, we need to redefine the concept of (/-classical OPS. Moreover, because the 
measures are supported on the real line, we will consider Q as the interval 1 [a, b]. 

DEFINITION 4.1. The sequence (p n ), where p n is a polynomial on x(s) of degree n for all 
n G No, is said to be a (/-classical OPS on the lattice x(s) if it satisfies the property of orthogonality 
(7) where 

(i) p(s) is a solution of the q-Pearson equation 

A(a( S )p(s))=T(s)p(s)Vx 1 (s). (15) 
1 If |6| < \a\ = oo (rcsp. |o| < \b\ = oo) then we write (a, b] (resp. [a, b)); if |6| = |o| = oo then we write (a, 6). 
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(ii) a(s) + ^t(s)Vxi(s) is a polynomial on x(s) of degree, at most, 2. 
(Hi) t{s) is a polynomial on x{s) of degree 1. 
(a(s), p(s), x(s)) is said to be (/-classical if it satisfies (i)-(iii). 

Notice that the above definition coincides with Definition 1.1 given in [4] for the g-linear lattices 
and generalizes it to the g-quadratic lattice. 

LEMMA 4.1. [24, (3.2.5), p. 62] For every polynomial ir and every integer k, ir(x k (s))+Tr(x k (s— 
1)) is a polynomial on Xk-i(s) of degree deg(7r). 

LEMMA 4.2. If (a(s), p(s),x(s)) is q-classical, then for every integer k and every polynomial 
ir, the function 

&l(p k {s),X k (s))(ir(x k+1 (s))) 

is a polynomial on x k (s) of degree deg(7r) + 1. 

Proof: By hypothesis, we know that a(x(s)) := a(s) + ^r(s)Vxi(s) and r(s) are polynomials 
on x(s) of degree, at most 2, and 1, respectively. From a direct calculation we get that 

. , . , u/1 , a(s + k) + t(s + k)Vxi(s + k) - a(s) 
^( Pk (s),x k (s))(l) = -i ' 

is a polynomial on x k (s) of degree l. 

Furthermore, for every integer k and every polynomial ir, we get 

sMs^m*^*))) = {ais+k)+ ' is+ " )Vxiis+k) K (x k+1 ( s) ) - ^ 



Vx k+1 (s) Vx k+1 (s) 

xir(x k+1 (s - 1)) = **( s + fc ) ' 7r(a; fe+ i(s)) - - ^ ir{x k+1 (s - 1)), 
\/x k+ i{s) \/x k+ i(s) 

being Q(s) := a(s) + r(s)Vxi(s). 

So, &i(p k (s),x k (s))(ir(x k+ i(s))) is a polynomial on x k (s) if and only if 

W^ 1(S), -^7W^- 1(S,) (16) 

is a polynomial on x(s). But we can rewrite this expression in the form 
0(s + |) -a(s- |) . . ., ct(s-|)/ / 

or, equivalently, 

9 '' + a , - >) '''- iw »^(' w '»'-' ( '- iw ')- 

Taking the arithmetic mean of the above expressions, using Lemma 4.1 as well as the relation 

@(s + |) - a(s - |) = gi(x(s))Vxi(s), 

where gi is a polynomial of degree 1, we deduce that (16) is a polynomial on x(s) of degree, at 
most, deg(-7r) + 1. Moreover, from a straightforward calculation the coefficient of x m+1 is 

a" 

ct q {m + 2k) t' + [m + 2k] q — / 0, m = deg(-7r), 

where a(x) = ^-x 2 + a'(0)x + a(0) and t(x) = r'x + r(0). Hence the result follows. □ 
Taking into account all these results, we can state the preliminary results related to our main 
Theorem. 



THEOREM 4.1. Let (p n ) be a q-classical OPS with respect to p(s) on x(s) such that 
x k (a)x l _ 1 (a)a(a)p(a) = x k {b)x l _ l {b)a{b)p{b) = 0, k,l > 0. 

Then, (A^p n+ i) is a q-classical OPS with respect to the function pi(s) on x±(s). 
Furthermore, if the condition (17) holds then the converse is also true. 



(17) 



REMARK 4.1. Observe that if, for instance, \a\ = oo (the case \b\ = oo is analogous) the 
relations regarding to a into (17) need to be replaced by 

lim x k (t)x l _ 1 (t)a(t)p(t) = 0, k,l > 0. 

Proof: Let (p n ) be a g-classical OPS with respect to p(s) on x(s). From Lemma 4.2 the function 
Qk(x(s)) := &i(p(s), x(s))(pk-i(xi(s))) is a polynomial on x(s) of degree k. Therefore for any 
n > k > 1 

6-1 

o= ^2p n {s)Q k {s)p{s)Vxi(s) 

s=a 
6-1 

= ^P„(s)(^i(p(s),x(s))(p fe _i(xi(s)))) p(s)V Xl (s) 

s=a 
6-1 

= J^Pn^V^pi^Pk^Xiis))^. 

s=a 

Taking into account the Leibniz rule, i.e. 

V(f(s)g(s)) = (Vf(s))g(s) + f(s - l)(Vg(s)), 
the following formula holds 



6-1 



6-1 



£ 9(s)(Vf(s)) = f(s)g(s + 1) - £ /(-)(A»( a )) 

^ — ' s=a— 1 ^ — ' 

s=a s=a 

Thus we get 

s=6— 1 

= Pn(s + l)/9l(s)p fe _l(xi(s)) - > (Ap n (s))pi(s)p fc _i(xi(s)) 

s=a— 1 ' — ' 
s=a 

k b-1 
s=b . ^ 

= p„(s)p(s)cj(s)p fc _i(x_i(s)) - \ (Ap„(s))pi(s)p fe _i(a;i(s)) (Aar(s) = Vx 2 (s)) 

s=a ' — ' 

s=a 

s=6 6-1 

= Pn(s)o-(s)p(s)pfe-i(x_i(s)) - V (A (1) p n (s))p fe _i(a;i(s)) pi(s)Vx 2 (s). 



(18) 



Hence (A^pn+i) is an OPS with respect to pi(s) on the lattice xi(s). 

Now, we prove that pi(s) satisfies a q- Pearson equation as (15) on the lattice x\(s). Indeed since 
(p n ) is a classical OPS we get 



A(a(s) Pl (s)) a(s + l) Pl (s + l) 



a(s) 



Pi(«) Pi( s ) 

= a(s + 1) + t(x(s + l))Vxi(s + 1) - cr(s) 

= a(s + 1) - 5(a) + -t(s + l)Vxi(s + 1) + -T(a)Vxi(s), 

where (j(s) = cr(s) + ^r(x(s))Vxi(s) is a polynomial on x(s) of degree at most 2, and after 
a straightforward calculation we deduce that this expression is equal to Ti(xi(s))Vx2(s) with 
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deg(fi) < 1. 

Moreover, using the last expressions we get 

<r(s) + \n{xi{s))Vx 2 {s) = - (a(s + 1) + a(s) + K(s + l)Vxi(s + 1) - -T(s)Vm(s)) , 
is a polynomial on X\(s) of degree, at most, 2 and hence the result holds. 

For the converse, we know that there exist two polynomials on xi(s), &(s) and n, of degree at 
most 2 and 1, respectively, such that 

A[a(s)p 1 (s)] = T 1 (x 1 (s))p 1 (s)Vx 2 {s), 



here a(s) := a(s) — ^T\(x\(s))Vx2{s). 
So, we only need to check that 



Vcj(s) 
Vx(s) 



and <t(s) + |r(s)Vxi(s) are polynomials on x(s) of degree at most 1 and 2, respectively, which 
is a direct calculation and hence the result follows. □ 

REMARK 4.2. Notice that in the case a = — oo (the case b = oo is analogous) one should write 
the formula (18) as 

]T g ( s )V(f( S )) = f(b-l)g(b)-(]imJ(t)g(t + l)\- ^ f(s)A(g(s)). 

s=— oo ^ ' s=— oo 

THEOREM 4.2. Let (p n ) be a q-classical OPS with respect to p(s) on x(s) such that 

x k (a)x l _ 1 (a)p(a) = x k {b)x l _ l {b)p{b) = 0, V k, I = 0, 1, . . . (19) 

If P-i = 0, then the sequence (<^ n (p_i(s), x_i(s))(l)) is a q-classical OPS with respect to the 
weight function P-i(s) = p(s — 1)/<t(s) on x_i(s). Furthermore, if the condition (17) holds then 
the converse is also true. 

The proof follows the same steps as in Theorem 4.1 taking into account the properties of A( m ) 
and the basic relations of p m (s) and x m (s) for every integer m. 

REMARK 4.3. 

(i) The relation between the statements of Theorems 4-1 and 4-2 follows from 

^i(p_i(s),x_i(s))o^ n (p(s),x(s)) =^ n+ i(p_i(s),x_i(s)), n = 0,1,2,... 

(ii) If the property of orthogonality 

b-l 

s=a 

and the boundary conditions (17) holds, then for every integer k we get 
b-k-l 



s=a~k 

where if k > 



J] (A^ Pn (x k (s))) (A^ Pm (x k (s)))p k (s)Vx k+1 (s) = 0, 
A(- fe ) :=<% k {p_ k {s),x- k {s)). 
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Now we can state the first main result of this paper. 

THEOREM 4.3. Let (p n ) be an OPS with respect to p(s) on the lattice x(s) defined in (14) and 
let a(s) be such that (19) holds. Then the following statements are equivalent. 

1. (p n ) is q-classical. 

2. The sequence (A^p n ) is an OPS with respect to the weight function p\(s) = a(s+ l)p(s + 1) 
where p satisfies (15). 

3. For every integer k, the sequence (& n (pk(s), Xfc(s))(l)) is an OPS with respect to the weight 
function Pk(s) where po(s) = p(s), Pk{s) = Pk-i(s + l)cr(s + 1), and p satisfies (15). 

4- (Second order linear difference equation): (p n ) satisfies the following second order linear 
difference equation of hypergeometric type 

A Vp n {s) ,&Pn{s) . , , , on s 
V*(s) + T ~Kx(s) + XnPn{3) = °' (20) 

where a(s) = a(s) + ^r(s)Vxi(s) and r(s) are polynomials on x(s) of degree at most 2 and 
1, respectively, and \ n is a constant. 

5. {p n ) can be expressed in terms of the Rodrigues Operator as follows 

p n (s) = BA(p(s),x(s))(l) = V V , . . (p n(s )) , (21) 

P(S) Vxl(S) \/X2{S) Vx n {s) 

where B n is a non zero constant. 

6. (Second structure relation,): There exist three sequences of complex numbers, (e n ), (/„), and 
(g n ), such that the following relation holds for every n > 0, with the convention p-\ = 0, 

, ,» Ap n+1 (s) A Pn (s) Ap n _!(s) 

^Pn(x(s)) = e n ^^ + f n ^y + <?n^p 
where ^# is the forward arithmetic mean operator: 

s /fU \. /(* + !) + /(*) 
^f(s) ■= - , 

e n 7^ 0, g n 7^ 7„ /or a// n > 0, and 7„ is the corresponding coefficient of the three-term 
recurrence relation [24] 

x{s)p n (s) = a n p n+1 (s) + n Pn{s) + TnPn-l(s), n > 0. 

REMARK 4.4. If u;e consider a q-linear lattice, i.e. either c\ = or C2 = 0, this result is 
"equivalent" to Theorem 1.3 stated in [4], because r(s)Vxi(s) is a polynomial on x(s) of degree 
2 and Xk(s) = q ak x(s) + 5k, where ctk, 5k are independent of s for every integer k. 

Assuming the theorem proved let us state and proof the second main result. 

THEOREM 4.4. Under the hypothesis of Theorem 4-3 the following statements are equiva- 
lent. 

(i) (Pn) is q-classical. 

(ii) (A«p n +i) is a OPS. 
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Proof: Of course (i)=>(ii) is a consequence of Theorem 4.1. 

(ii)=>(i): This proof is analogous to the proof given by W. Hahn in [15]. We know that the SPO 
(p n ) satisfies the TTRR 

Pn(x(s)) = (x(s) - a n )p n -i(x(s)) - b n p n - 2 (x(s)) , n > 1, 

with initial conditions p-\ = and po = 1. Thus applying the (forward) arithmetic mean process 
we get 

A^p n (x(s)) = JZp n - X {x{s)) - b n A^p n _ 2 (x(s)) + (a q (2) Xl (s) - a n )A^p n ^{x(s)), (22) 

with a n = a n - c 3 (q)(l - a q {2)). 

Moreover, by hypothesis, the sequence of monic polynomials {[n]~ l A^p n ) also satisfies a 
TTRR, i.e. 

J-A( 1 )p n (x( S )) = (x 1 ( S )-a^) ?; ^A(V-i(x( S ))- T -^ r A(V~2(x(,)), n > 1. (23) 

Combining (22) and (23) to eliminate A^p n _ 2 (resp. A^p n ) we get 

{A ~ k) ^Ms)) = ((fgK - j^h(s) ~ ^ + £t) A(Vi(xW) 

(24) 



and 

= (S^fe^) - A- q + a) A( V-i(*( s )) 

- (A - ft) A(V- 2 (x( S )). 
Setting n + 1 instead of n in the above expression the following expression fulfills: 



(25) 



In order to simplify further calculations in this proof we denote by e n := ]^ 2 ] |nT> ^( x > n ) = 

[n+\] q -[n-l] q _ <+i , Qn+i j l \ _ / Qg(2)6^ _ fe n \ _ , fena^ 

2[ra,][n+l], X M + [n+l] 9 allCL °y H^.^J - ^ [„_ 2 ] q [n-l]J X [n-2], "+" [n-1], • 

Taking into account (24) and (25) in order to eliminate A^p n we obtain 
s b' 

^-T-^p n (x(s)) = 2 - r A:(a;i(s);n)^p„_i(a;(s)) + fc(a;i(s);n)(Z(a;i(s);n)-e„ + ie n )A (1) p rl _i(a;(s)). 



[n+l] g [n-2] g 

(26) 



By using (25) and (26) to eliminate p n we get 

,#(l(i-i(s);n)Vp„-i(x(s 

») , & .*((VgOOMfrn-l(s(s-l))) 

[n+l] g 6; V(fc(xi(s);n)^p„_i(x(s))) . , V((/(a;i(s);n)fc(a; 1 (s);n)-e„e n+ i)A( 1 )p n _i(a;(s))) 
- [n-2j, WTTs) h [ n+1 ]« VST^l • 

On the other hand, combining (24) and (25) to eliminate A^p n -i we get 

b j^^p n _ 1 (x(s)) = l -^>l^p n (x(s)) - (A;(x 1 ( S );n)/(x 1 ( S );n) - e n e n+1 )A^p n (x(s)). 

(28) 
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By using the expressions (25) and (28) to eliminate p n -\ we get 

1 J?{{Vx(s))Ji(pn{x{s-l))) _ ^(k(x-i(s);n)Vp n (x(s))) 

[n-2] q V(f(ai(s);n).^p ra (a:(s))) _ [n-2] q V((K^i(g);»)fc(a:i(^);")-ene T1+ i)A( 1 )p n (a;(s))) 
- V n [n+l\ q WxT(s) b' n V^TIS) • 

The following result has a technical role and it will be used in to proof of this Theorem. 
LEMMA 4.3. For any two polynomials tt and the following relation holds: 

MM^^Msm = vfrfrwy (,))) + po.) <Xl{8)mx{8)) 

- WxlJS) M Vx( s ) + v{n{xi{s))l\x{s)) Vxi[s) — . 

Moreover, the functions ■^( n ( x ^^))^ 7x ( s )) andV(Tr(xi(s))Ax(s)) are polynomials onx(s). 



The proof of this result is straightforward and we leave it as an exercise for the reader. 
Therefore, using Lemma 4.3 and doing some extra calculations, we obtain that the expressions 
(27) and (29) become the following two second order homogeneous linear difference equations: 

2 (x( g ),n)v( 1 )A(V(^(^)+0i(x( g ),n)^ V ^ ( ) S)) + fo(nK(s(s)) = 0, 

<j> 2 (x(s), n)V^A^ Pn ^(x(s)) + <Mx(s), n)^ VPn ~^ s)) + 4> (n)Pn-i(x(s)) = 0, 

where the indices indicate the degree of the polynomial coefficients. 

Moreover, replacing n by n + 1 in the last expression we obtain the second order linear 
difference equation: 

<h{x(s),n + l)V^A^ Pn (x(s)) + Mx(s),n + l)Jt Vp ™^ ] + + l) P n(x(s)) = 0. (30) 

So, we have two difference equations for p n which only differ in a constant factor. 
Therefore p n satisfies the following two second order linear difference equations: 

V2 (x(s))V^A^ Pn (x(s)) + yi (x( g ))^ V ^ (s)) + Mx(s))p n (x(s)) = 0, 

Mx(s))V {1) A^p n (x(s)) +^ 1 (x( S ))^ V ^y ) +Mx(s))p n (x(s)) = 0, 

where, again, the indices indicate the degree of the polynomial coefficients. Denoting by a(s) := 
(p2{x(s)) + ipi(x(s))Vxi(s), ((s) := ip2(x(s)) + ipi(x(s))Vxi(s), the above expressions can be 
rewritten as 

a(s)V^A^p n (x(s)) + ^x{s))A^p n (x(s)) + m (n) Pn (x(s)) = 0, 

C(s)V«A( V(*(s)) + Mx(s))A^ Pn (x(s)) + Mn)Pn(x{s)) = 0. 

Notice that a and £ are not always polynomials on x(s). 
So, 

(^i(x( S ))C(s) - a(s)Mx(s)))A( 1 ) Pn (x( S )) 
+(<Po(n)((s) - a(s)ip (n))p n (x(s)) = 0. 

If n > 3, we get that p n and A^p n , has a common zero but p n has not no multiple zeros. Hence, 
after dividing by the common factor, the polynomial coefficients of the difference equations for 
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Pn+i, Pn, and p n -i are the same. Taking this into account follow the coefficients do not depend 
on n, so p n (x) satisfies the second order linear difference equation 

<P(x(s))V^AWp n (x(s)) + V>(x( S ))-^ V ^ (s)) + KPn(x(s)) = 0, 

where degcj) < 2, degtft < 1, and A n = — [n] 9 (a 9 (n — 1)6^, + [n — l]a^) depends on n since p„ is 
a polynomial, being and fry, the leading coefficients of <^> and Y>, respectively. 

Hence, according to the 4th condition of Theorem 4.3, (p n ) is (/-classical. □ 
Let us prove now Theorem 4.3. From Theorem 4.1, Theorem 4.2, and its Corollary, as well 
as Remark 4.3 we know that (1) -> (2) -> (3) -> (1). 

PROPOSITION 4.1. ((l)-> (4)) // (p n ) is a q-classical OPS with respect to p(s) satisfying 
the boundary relations (19) then (p n ) satisfies the second order linear difference equation of 
hypergeometric type (20). 

Proof: By Lemma 4.2 the function &i(p(s), x(s))(A^p n (x(s)) is a polynomial on x (s) of degree 
n. On the other hand, by Theorem 4.1 (A^p n ) is a (/-classical OPS with respect to pi(s) on 
xi(s). Thus, for < k < n, 

6-1 

Y,Qk(x(s))^ 1 (p(s),x( S ))(AWp n (x( S )))p( S )Vx 1 ( S ) 

s=a 

= b '£Q k ( X (s))V(p 1 (s)AWp n (x(s))) 

s=a 

t. b-1 
s=b 



(18) 



Q k (x(s))a(s)p(s)A^p n (x(s-l)) ' A^(p n (x(s)))A^(Q k (x(s))) p 1 (s)Vx 2 (s) 

Q CI ' * 



By hypothesis and Theorem 4.1 the last expression vanishes and, as a consequence, there exists 
a non zero constant A n , independent of s, such that 

M l {p{s),x{s)){A^p n {x{s))) = -X nPn (x(s)), 

since (p n ) is an OPS with respect to p(s) on x(s). 

Finally, using that the above expression is indeed Jif q (see (13)) it is clear that the above 
expression (20) and hence the result holds. □ 

(4) — >(5): This is a well-known property. See, for instance, §3.2.2 in page 64 of [24]. 

(5) — >(1): By setting n = 1 in the Rodrigues formula we obtain the Pearson type equation (15) 
that the weight function p(s) satisfies. 

PROPOSITION 4.2. ((6) <-> (1)). Let (p n ) be an OPS with respect to p{s) on the lattice x(s) 
defined in (14) and let a(s) be such that (19) holds. The sequence (p n ) is q-classical if and only 
if there exist sequences of complex numbers, (e n ), (f n ), and (<?„), such that the following relation 
holds 

J/ v ( x ( s )) - e A Pn+l(x(s)) f &Pn(x(s)) A Pn ^(x( S )) 

M Vn {x{s))-e n A ^ +U A ^ +g n , 
with the convention p-\ = 0, where e n ^ ; g n ^ 7 n for all n > 0. 

Proof: If (p n ) is an OPS then it satisfies a TTRR, i.e. there exist three sequences of complex 
numbers, (a n ), ((3 n ), and (7 n ), with j n ^ 0, such that 

x(s)p n (x(s)) = a n Pn+i(x(s)) + n p n {x(s)) + j n Pn-i (x(s) ) , n > 0. (31) 

If (p n ) is g-classical, then (A^p n+ i) is an OPS, and therefore there exist sequences of complex 
numbers, (ai 1 '), {Pn^), and (7n^), with / 0, such that 

Xl {s)A^p n {x{s)) = a ^AW Pn+1 (x(s)) + tipAW Pn (x(s)) + 7 «A( 1 V n _ 1 (x( S )). (32) 
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But 

A^(x(s) Pn (x(s))) = (^x(s))A^ Pn (x(s)) + J(p n (x{s)). (33) 
So combining (31), (32), and (33) we get 

JZ Pn {x{s))=(a n - [ ^aWy^p n+1 (x(s))+p n AW Pn (x(s)) + ( 7n - ^^A^-i 



where 



Moreover, the coefficient of A^p n -i is different from 7n because 7 n^ 7^ 0. 

Conversely, if there exist sequences of complex numbers, (e n ), (/„), and (g n ), such that the 

following relation holds 

assuming p_i = 0, then from (33) we get 

^ Xl (s)A^p n (x(s)) = (a n - e„)A( 1 Wi(^(s)) + A( 1 )p n (x( S )) + ( 7n - ^A^-i 



where 



^ =(3n-C 3 -f n + l -fc 3 . 



By hypothesis # n / 7n , hence by Favard's Theorem (A^^ n+ i) is an OPS, and by Hahn's 
Theorem the result holds. □ 



5. The examples 

5.1. The Askey- Wilson Polynomials 

The Askey- Wilson polynomials, which were introduced by R. Askey and J. Wilson, are located in 
the top of the q- Askey tableau [16]. These polynomials can be written as a basic hypergeometric 
series 



p n (x(s);a,b,c, d\q) 



(ab; q) n (ac; q) n (ad;> 



■4^3 



~n—l 



q n , abcdq n aq" , aq 
ab, ac, ad 



where a, b, c, d, ab, ac, ad, be, bd, cd {q m : m £ Z}. 



REMARK 5.1. Notice that although this family does not satisfy a property of orthogonality 
(7) it is orthogonal with respect the linear functional u AW which has the following integral 
representation [16, (3.1.2), p. 63] 



(u AW ,p)= ^ p(x)u(x)dx 



p(x) 



1 v 7 ! — x 2 



(e 2W ;q) c 



(ae i9 ,be ie ,ce ie ,de ie ;q) c 



dx, x = cos 9. 



Furthermore, such a functional fulfills the distributional equation 

AW(a AW u AW ) = Pl (x(s); a, b, c, d\q)u AW , 

and it is straightforward to check that (a AW (s),u>(s),x(s)) is q-classical, where x(s) = ^(q s + q^ s ), 

i.e. c\ = C2 = \ and c 3 = 0, a AW (s) = — K 2 q q~ 2s+ z (q s — a)(q s — b)(q s — c)(q s — d), and hence we 
can apply Theorem 4-3 to Askey-Wilson polynomials. 
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These polynomials are the polynomial eigenfunctions of the linear difference operator [16, (3.1.6)] 
jg^w = ^_ Law { _ s) A _ aAW V \ 

Notice that if r AW (x(s)) = pi(x(s); a, b, c, d\q) a straightforward calculation yields to 

a AW (-s) = a AW (s) + r iw (x(s))Vn( S ). (34) 

Here, as before, K g = <p — q~2 . With these parameters, let us consider the following function 

p AW (s) = q~ 2s2 (a,b, c, d; q) s (a,b, c, d; q)- s , 
where (a; q)o = 1, (a; q)k = (1 — a)(l — aq) ■■■(! — ag fc_1 ), k > 1, and 

9) - fc = (1 - ag -i)(l - ag -2) ••• (1 - ag -*) ' fe " 1 ' 
LEMMA 5.1. XTie function p AW (s) satisfies the Pearson-type equation 

a AW (s + l)p(s + l)=a AW (-s)p(s). (35) 

Proof: 

p AW (s + 1) _ g- 2 ^ 1 ) 2 («, 5, c, d; g) a+1 (a, b, c, d; q)^ 
p AW (s) q- 2s2 (a,b,c,d;q) s (a,b,c,d;q)- s 

q -4s-2^ _ Qqs^ _ _ _ dq s^ 



(1 - ag- s - 1 )(l - 6g- s - 1 )(l - cq- 8 ' 1 )^ - dq- 3 - 1 ) 

q 2s {q~ s - a){q~ s - b)(q- s - c){q- s - d) a AW (-s) 



- q -2s-2( q s+l _ a )( ?s +i _ b )( g .+i _ c ^ q s+i _ d ) a AW( s + !) 

□ 

Notice that from (34), the equation (35) becomes the Pearson-type equation (15). 
Taking into account that p n (s) = p n -i(s + l)o"( s + !),"-> 1, with po(s) = p(s), a straight- 
forward calculation gives 

Pn W (s) = 4 n q- 2s2 ~ 2sn ~ n2+3 * n (a,b,c,d-,q) s+n (a,b,c,d;q)- s . 
Then the Askey- Wilson polynomials can be written for every nonnegative integer n as 

/ / s , us B n nl n q 2s2 (n) (a,b,c,d;q) s+n (a,b,c,d;q)- s 

p n (x(s):a,b,c,d\q) = -. — ; ; — r-^ — ; ; — r — V v ; — - , 

y V ' W (a,b,c,d;q) s (a,b,c,d;q)_ s q 2s*+2sn+n*-ln 

_ _ n(3n-5) 

where B n = 2 n K„ n q * 

Furthermore, since this family fulfills the following difference relation [16, (3.1.8)] 

ls! A) p n (x{s);a, b,c,d\q) = 2[n] q (l - abcdq n ' 1 )p n -i(xi(s); aq^ ,bq^ , cq^ , dq^\q), 

and the coefficients of the second structure relation are 

2[n],(l - abcdq n - 1 ) 2 - [2] q [n + 1],(1 - abcdq 71 ) 2 
6n ~ 4[n] 9 (l - abcdq n - l ){l - abcdq 2n - r )(l - abcdq 2n ) ' 

fn = ~^( a ~ a ~ 1( T l ) ~ 2 (^ A n(a,b,c,d\q) + C n (a,b,c,d\q) 

-^(An-i (aq^ , bq^ , cq^ , dq^ \q) + C n _i (aq^ , bq^ , cq^ ,dq*\q))j , 
(1 - abcdq^- 2 ^){l - abcdq 2n ~ 2 )(l - abcdq 2 ^ 1 ) ( 2[n] q (A n _ 1 C n )(a,b,c,d\q) 



9n 



4[n] g (l - abcdq 11 - 1 ) 2 \ (1 - abcdq^- 2 ^ 

[2]g[n - l] g (^[ n -2] g C l »-i)(ag^,^g^,cg^,dg^|g ) 
(1 - abcdq n ~ 1 ) 2 
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REMARK 5.2. Notice that in [19] T.H. Koornwinder has obtained a different explicit structure 
relation for Askey-Wilson polynomials by using the difference operator theory. 



5.2. The (/-Racah polynomials 

We consider the g-Racah polynomials u^{x{s),a,b) q on the lattice 

x{s) = [s] q [s + l] q = q^K~ 2 q s + q-^K q 2 q~ s - [2] g ^ 2 . 



They were introduced in [20, 24] and deeply studied in [6]. These polynomials can be written in 
terms of a basic hypergeometric series [6] as follows 



/ Q—n „a+/3+n+l „a-s a+s+1 
<' f '(x(s),a,b) q =D nA <p 3 [ „/3+l n a+b+a+l 



Q,q] , 



where 

q -^(2a+a+t3+n+l)^ q a~b+l. q ^ n ^+l. q ^ q a+b+a+l. q ^ n 



D n = 



n q ]n (q;q)n 



Observe that from the above formulas the polynomials u^{x{s) } a, b) q are multiples of the 



standard q-Racah polynomials £itf n (fj,(q b+s );q a ,q l3 ,q a ~ b ,q~ a ~ b \q). 

These polynomials are the polynomial eigenfunctions of the second order homogeneous linear 
difference hypergeometric operator 

jm R = —1— (a qR (-s - 1)—=^- - a qR (s)-^--) , 
q Vxi(s) V J Ax(s) w Vx(s)y' 

where 

o- qR (s) = -^L^ {q s_ q a ){q s_ q -b ){q s_ q p-a ){qS _ qb+ ^ = [ s - a ] q [s + b] q [s + a - (3} q [b + a - s] q . 
K\q 2 

Furthermore, taking into account the ^-difference relation 

A^u^(x(s), a, b) q = [a + (3 + n + l^'^iW.a + \,b- \) q . 
The coefficients for the second structure relation are 



2[n + a + (3 + l]g - [2] q [n + a + (3 + 2} q 
2[2n + a + j3 + l] q [2n + a + (3 + 2] q [H J<? ' 



/„ = (3 n (a, b, a, (3) - [ -^(3 n (a + \, b - \, a + 1, f3 + 1) 



[2Js 
2 

[a + b + a + n] q [a + b — (3 — n] q [a + n] q [(3 + n] q [b - a + a + (3 + n] q [b - a - n] q 
9n ~ 2[a + + 2n]Ja + /? + 2n + l] q [n + a + /? + 1], 

x(2[n + a + /3+l] g - [2] ? [n + a + f3] q ). 

5.3. The (/-Meixner polynomials 

This family of g-polynomials has the following representation as a basic hypergeometric series 



(bq;q) n (-c) n ( q~ n ,x 



M n (x; b, c; q) = - — - 2 2 <fi 



bq 



1 n+1 \ , ^ - 
, x(s)=q =x. 
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They are the polynomial eigenfunctions of the linear difference operator of hyper geometric type 
[16, (3.13.5)] 



— — ( (x — l)(x + be)- — ;-r — q 1 c(x — bq)—^-—\ , 
i(s) V Ax(s) H 1 4; Vx(s)y 



9 Vx 

In this case 

_ [2]Jn + l] _i 
2[n], 9 ' 

1 _ n n -2n-l rr 2n 
f n = —1+ q —— ( c - ebq n+1 ){2 - q - q *) + l—(c+q n )(2-q- q 2 + q n - q^ 1 ), 

9n = Cg " 4W+1(1 9 7^ )(C + gW) (2N 9 (l " 1 n ) - (<Z + 9 2 )[n - 1U1 - q^ 1 )). 

z l n \q 

5.4. The Al-Salam & Carlitz polynomials I and II 

The Al-Salam & Carlitz polynomials I (and II) appear in certain models of (/-harmonic oscillator, 
see e.g. [8, 9, 10, 23]. They are polynomials orthogonal on the g-linear lattice x(s) = q s = x, 
defined [16] by 

,-n „-l qx 

2<Pl ' 



Ui a \x-q) = {-aTq^) W A q 



These polynomials are the polynomial eigenfunctions of the linear difference operator of hyper- 
geometric type [16, (3.24.5)] 

/ A - (x - 1)(x " a) v^))' 



ACI - 



a- 



q Vxi(s) V Ai(s) 

In this case [16, (3.24.7)] A^lli a \x;q) = q^~[n] q U^} 1 (x;q), so the second structure relation 
coefficients are 



V.- 1 T H V" ' ^ q f /-, . 

e n = 1 z-p; -, fn= (1 + a 

2q[n] q 



2[n], 



5 



5n = -^j- (2[n]g3( g » - 1) - [ 2 ]Jn - 1],^ - 1 

Taking into account that [16, p. 115] 

V^(x;q) = Ui a \x;q- 1 ), 

all the information related to the Al-Salam & Carlitz polynomials II can be deduced from the 
information for the Al-Salam & Carlitz polynomials I in a simple way by changing q to q^ 1 . 
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